clc;
clear;
A=[16 2 3 13;
    5 11 10 8;
    9 7 6 12;
    4 14 15 1];
[L1,U1]=lu(A)
disp('-----------------------------')
[L,U,P]=lu(A)
[L2,U2]=lu(sym(A))
